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Abstract 



Physical Review E Vol. 58, No. 6, p. 7013-19 



The noise of signals or currents consisting from a sequence of pulses, ele- 
mentary events or moving discrete objects (particles) is analyzed. A simple 
analytically solvable model is investigated in detail both analytically and nu- 
merically. It is shown that 1/f noise may result from the statistics of the pulses 
transit times with random increments of the time intervals between the pulses. 
The model also serves as a basis for revealing parameter dependences of 1/f 
noise and allows one to make some generalizations. As a result the intensity 
of 1/f noise is expressed through the distribution and characteristic functions 
of the time intervals between the subsequent transit times of the pulses. The 
conclusion that 1/f noise may result from the clustering of the signal pulses, 
elementary events or particles can be drawn from the analysis of the model 
systems. 

PACS number(s): 05.40.+j, 02.50.-r, 72.70.+m 
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I. INTRODUCTION 



The omnipresence of 1/f noise is one of the oldest puzzle of the contemporary physics. 
During more than 70 years since the first observation by Johnson, long-memory processes 
with long-term correlations have been observed in many types of systems from physics, 
technology, biology, astrophysics, geophysics and sociology (see and references herein). 
Recently 1/f noise was discovered in human cognition [f|, human coordination || and even 
in distribution of prime numbers H]. 

Despite of the widespread occurrence of fluctuations of signals and variables exhibiting 
1/f 5 (<5 ~ 1) behavior of the power spectral density S (/) at low frequencies in large di- 
versity of systems no generally recognized explanation of the ubiquity of 1/f noise is still 
proposed. Physical models of 1/f noise in some physical systems are usually very specialized, 
complicated (see and references herein) and they do not explain the omnipresence of 
the processes with 1/ f s spectrum PHTH- Note also some mathematical analysis [IT ], models 



and algorithms of generation of the processes with 1/f noise fT2}-|T4||. These models also ex- 
pose some shortcomings: they are very specific, formal (like "fractional Brownian motion" 
or half-integral of a white noise signal) or unphysical. They can not, as a rule, be solved an- 
alytically and they do not reveal the origin as well as the necessary and sufficient conditions 
for the appearance of 1/f type fluctuations. 

In such a situation the simple analytically solvable model system generating 1/f noise 
may essentially influence in the reveal of the origin and essence of the effect. Here we 
present a model which generates 1/f noise in any desirably wide range of frequency. Our 
model is a result of the search of necessary and sufficient conditions for the appearance of 
1/f fluctuations in simple systems affected by the random external perturbations initiated 
in |lj| and originated from the observation of the transition from chaotic to nonchaotic 



behavior in the ensemble of randomly driven systems [16 . Contrary to the McWhorter 



model [DJ based on the superposition of large number of Lorentzian spectra and requiring 
a very wide distribution of relaxation times, our model contains only one relaxation rate 7 
and can have an exact 1/f spectrum in any desirably wide range of frequency. The model 
may be used as a basis for the checking of assumptions made in the derivation of 1/f noise 
spectrum for different systems. Furthermore, it allows us to make a heuristic presumption 
for the generalizations of the theory of 1/f noise. Numerical simulations and comparisons 
with analytical results confirm this supposition. 



II. THE MODEL 

In many cases, the intensity of some signal or current can be represented by a sequence 
of random (however, as a rule, mutually correlated) pulses or elementary events (t — t^). 
Here the function Ak (ip) represents the shape of the k pulse having an influence to the signal 
I (t) in the region of transit time t^. The signal or intensity of the current of particles in 
some space cross-section may, therefore, be expressed as 

/(*) (i) 

k 

It is easy to show that the shapes of the pulses influent mainly on the high frequency, / > At p 
with At p being the characteristic pulse length, power spectral density while fluctuations of 



2 



the pulse amplitudes result, as a rule, in the white or Lorentzian but not 1/f noise [18 . 
Therefore, we restrict our analysis to the noise due to the correlations between the transit 
times tfc. In such an approach we can replace the function A k (t — tk) by the Dirac delta 
function 5 {t — t k ) and the signal express as 

/(*)=£*(*-**)■ (2) 

k 

This model also corresponds to the flow of identical objects: electrons, photons, cars and so 
on. On the other hand, fluctuations of the amplitudes A k may result in the additional noise 
but can not reduce 1/f noise we are looking for. 
The power spectral density of the current (2) is 



S(f) 




max 



k 

k—k ■ 



/ e 




^max k 



where T is the whole observation time interval, k m \ n and /c max are minimal and maximal 
values of index k in the interval of observation and the brackets (...) denote the averaging 
over realizations of the process. 

In this approach the power spectral density of the signal depends on the statistics and 
correlations of the transit times t k only. It is well known that sequence of random, Poisson, 
transit times generates white (shot) noise [[Il|. The sequence of transit times t k with random 
increments, t k = t k _i + f + ae k (where f is the average time interval between pulses, 
{e k } denotes the sequence of uncorrelated normally distributed random variables with zero 
expectation and unit variance, i.e. the white noise source, and a is the standard deviation 
of white noise) results in the Lorentzian spectra [jnj . Here we will consider sequences of the 
transit times with random increments of the time intervals between pulses, r k = r k -i + <JSk, 
where Tk — tk — tk-i- It is natural to restrict in some way the infinite Brownian increase or 
decrease of the intervals Tk, e.g. by the introduction of the relaxation to the average period 
f rate 7. So, we have the recurrent equations for the transit times 

= tk-l + Tk, / 

= Tk-i - 7 {r k -.i -f)+ae k . 

The simplest physical interpretation of the model (4) corresponds to one particle moving 
in the closed contour with the period of the drift of the particle round the contour fluctuating 



(due to the external random perturbations) about the average value r [19 . 



III. SOLUTIONS 

An advantage of the model (4) is that it may be solved analytically. So, an iterative 
solution of Eqs. (4) results in an expression for the period 

r, = r + (r -r)(l- 7 f + ^(l- 7 ) H ^ (5) 
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where r is the initial period. The dispersion of the period is 



^(*) = (7*)-fo> 



0" 



l-(l- 7 ) 



2 A- 



27(1-7/2) 



cr 2 A;, 2£ry < 1 
<7 2 /2 7 , 2A; 7 > 1. 



(6) 



Therefore, after the characteristic transition to the stationary process time, tt r — f/"f, the 
dispersion of the period approaches the limiting value a 2 = a 2 /2 / ~f. 

After some algebra we can also obtain an explicit expression for the transit times tk 

(*>1), 



h — to + kr + 



1-7 

7 



l"(l-7) fc ] (T-o-r) + ^E[l-(l-7) 

7 i=i 



k+l-l 



(7) 



where t is the initial time. The dispersion of the transit time tk is 



°i(k) = (ti )-(tk) 



,2 <r_ 

r 



fc-2 



1-7 

7 



l-(l-7)" +(l-7) : 



l-(l-7) 
1-(1- 7 ) S 



2fe 



(8) 



= f a 2 (k/6 + k 2 /2 + k 3 /3 + ...), 2 7 £;<1, 
~ \ (a/7) 2 (A; - 3/2 7 + 5/4 ± ...) , 2 7 A; > 1. 

At /c ^> 7 _1 Eq. (7) generates the stationary time series. The difference of the transit times 
t k+q and tk in Eq. (3) for r = f or 27/c 1 is 



- tk = fq + a -\ [1 - (1 - 7f] E (1 " 7) fc+1 "^ + E I 1 " (! " 7) 



The dispersion of this times difference equals 



where 



l=k+l 



a 



(t k+q -t k ) ) - t q = — g (q) 



:i},q>0. 
(9) 

(10) 



9 (?) = ^ {[1 - (1 - 7)f E(l- 7) 2 ' + E [!-(!- 7) 



g > 0. 



1=1 



=1 



Summation in Eq. (11) results in 



9 (?) = -2 

7 2 



9- 



(1 



[l-(l-7) 4 



i-(i- 7 r 



{2 + [l-(l-7)1(l-7) 2fe+1 } 



At 7a < 1 

while for 27/c > 1 we have 



(2k + 1) g 2 + q/3 + 2q 3 /3, 2-fk < 1, 

( 7 + 1) <? 2 + i? - §? 3 , 27A; » 1, 



(11) 



(12) 



(12a) 



Note that for g < we should replace in Eq. (9)-(13) q by |g|and k by k — \q\. Therefore, 
the function g(q) at k — \q\ ^> 7 _1 is even, i.e. g{—q) = g(q)- 

Substituting Eq. (9) into Eq. (3) and replacing the summations in the exponents by 
the multiplications of the exponents we have the following expression for the power spectral 
density of the current 

S(f)= hm /|E ""if* e^nexp{^[l-(l- 7 n(l- 7 ) fe+1 -^,} 

T^OO \ k g=fc min -fc 1 = 1 J (g fl ) 

"ff exp{i^[l-(l- 7 )^ +1 -'] £ ,}\. 

The average over realizations of the process coincides with the average over the distri- 
bution of the random variables S\. Using the fact that random variables E\ are independent 
and mutually uncorrected we can fulfill the averaging over every random variable S\ inde- 
pendently. Therefore, Eq. (3a) may be rewritten in the form 

S (f) = hm | E ft (exp [1 - (1 - if] (1 - 7)^ 

T-oo k>q i =1 \ l 7 /gjx 

[l-(l-,r>}>. 

The result of the averaging of the exponent exp {ice;} (with c being a constant) over the 
normally distributed random variable E\ with zero expectation and unit variance is 

+00 



—00 



Therefore, after the averaging over the normal distribution of the random variables Si Eq. 
(3b) takes the form 

S (/) = hm I E ft exp [1 - (1 - 7 )f (1 - 

*ffexp{-^[l-(l- 7 r»«-'] 2 }. 

Z=fc+1 I J J 

Transition in Eq. (3c) from the multiplications of the exponents to the summations in 
the exponents and transformations in analogy with Eq. (11) of the two sums' summation 
indexes I — > k + 1 — I and I — > k + q + 1 — I, respectively, yield according to Eq. (11) the 
final expression for the power spectral density 

S (/) = hm ^ £ e *2*/^-*Wff(9). (14) 

k,q 
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Since the expansion of the function g (q) in powers of 7 \ q\ <C 1 at 2^k ^> 1 according to 
Eqs. (12) and (13) is 

g(q) = -q 2 -\\q\ 3 + lq 2 ±...., (15) 

for f <ti f T = (27rf) _1 and / < f 2 — 2^Py/na we can replace the summation in Eq. (14) by 
the integration 

+00 

S (f) = 21 J e i2nffq ^ 2f2a29 ^dq. (16) 
— 00 

where I = lim (k max — k min + 1) jT — f~ l is the averaged current. 

T^oo 

Furthermore, at / f± — '-f'l 2 jixo it is sufficient to take into account only the first term 
of expansion (15), g (q) = q 2 / r y. Integration in Eq. (16) hence yields to 1/f spectrum 



S(f) = 2I J exp 



-00 



,2 



l 2nfrq - { -^q 2 

7 



dq = I 2 ^, /i</</ 2 ,/r (17) 



where an is a dimensionless constant (the Hooge parameter) 



a H = ^=Ke~ K \ K=-^ = ^. (18) 

V2a T cr 



Using an expansion of the function g(q) at jq 1 according to expression (13), g(q) = 
2q/ r y 2 , we obtain from Eq. (16) the Lorentzian power spectrum density for / < f\ 

S (/) = 21^- , = I 2 4T f - . (19) 

K ' f 2 7 2 1 + (nfa 2 /f 7 2 ) 2 1 + r> 2 1 ' 

Here u = 2irf and r re / = D t = a 2 j2f^ 2 is the "diffusion" coefficient of the time tk according 
to Eqs. (7) and (8). The model is, therefore, free from the unphysical divergence of the 
spectrum at / — > 0; for / <C f = f^/na 2 = l/27ir re i we have from Eq. (19) the white 
noise 

S(f) = I 2 (2a 2 /f 7 2 ) . (20) 

Therefore, the model containing only one relaxation rate 7 for sufficiently small parameter 
7 can produce an exact 1/f- like spectrum in any desirably wide range of frequency, f\ < 
f < f'2, fr- Furthermore, due to the contribution to the transit times tk of the large number 
of the very separated in time random variables, E\ (I = 1,2, ...k), our model represents a 
'long-memory' random process. 



IV. GENERALIZATIONS AND NUMERICAL ANALYSIS 

Eqs. (16)-(20) describe quite well the power spectrum of the random process (4). As an 
illustrative example in Fig. 1 the numerically calculated power spectral density averaged over 
five realizations of the process (4) is compared with the analytical calculations according to 
Eqs. (16)-(20). The analytical results are in good agreement with the numerical simulations. 
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Note, that analytical results predict not only the slope and intensity of 1/f noise but the 
frequency range /i -j- f%, f T of 1/f noise and intensity of the very low frequency, / <S /o, 
white noise (20) as well. 

This model may also be generalized for the non-Gaussian and for the continuous per- 
turbations of the systems' parameters resulting in the fluctuations of the period r. So, for 
perturbations by the non-Gaussian sequence of random impacts {ek} with zero expectations 
Eqs. (1)— (13) remain valid. Only the result (14) of the averaging over realizations of the 
process in the case of the non-Gaussian perturbations may have different form. Consider 
now such a situation in more detail. 

The power spectral density (3) may be rewritten in the form 

S (f) = 21 (j2 e i2wfrkiq)q ^ (21) 
where the transit times tk+ q and tk difference is expressed as 

k+q 

t k+q -t k =J2 T k = T k(q)q, q>0 (22) 

l=k+l 

and the brackets denote the averaging over the time (index k) and over the realizations of the 
process. Here r k (q) = {tk+q — tk) /q is the averaged time interval between the subsequent 
transit times in the time interval tk-^tk+q- Note that for the slow (diffusive-like) fluctuations 
of the averaged interval T\. (q) with the change of the index k Eq. (22) is valid also when 
q < , i.e. t k+q -t k = T k - q (q) q ~ r k (q) q, q < 0. At 2nfT k (q) « 1 we may replace the 
summation in Eq. (21) by the integration and do not take into account the dependence of 
Tfc (q) on q. In such a case Eq. (21) yields 

; +00 , +00 

S(f) = 21 1 j e l27TfTkq dq\ = 21 J (e i2nfrkq ) dq. (23) 

oo ' — oo 

Here the averaging over k and over the realizations of the process coincides with the averaging 
over the distribution of the periods Tk, i.e. 

+oo 

(j**f<rn^ = J e i2wfqT ^J (r) dr = Xr (2irfq) (24) 

— oo 

where ip (r) is the periods distribution density and Xt ($) is the characteristic function of 
the distribution ip (r). Taking into account the property of the characteristic function 

+oo 

J X T (#)d# = 27T^(0) 
— oo 

we have from Eqs. (23) and (24) the final expression for the power spectral density 

S(/) = 2/^(0)//. (25) 



7 



Substituting into Eq. (25) the value ^ (0) = exp (— f 2 /2of) / ^/2na T for the Gaussian 
distribution of the periods tj, we recover the result (17)— (18). 

Since different processes result in the Gaussian distribution it is likely that perturbation 
by the non-Gausssian impacts {sk} in Eq. (4) yields nevertheless the Gaussian distribution 
of the periods r k . For the demonstration of such statement and validity of the approach (21)- 
(25) we have performed numerical analysis of the model (l)-(4) for different distributions of 
the perturbations {sk}- Figures (2) and (3) represent the calculated power spectral densities 
for the rectangular (uniform) and asymmetric xi distributions of the sequence {e^} with 
zero expectations and the same variances and other parameters like those for Fig. 1. We 
can notice only the slight dependence of the spectra on the distribution function of the 
perturbing impacts {e^} with the same expectations and variances. These results confirm 
also the presumptions made in the derivation (21)-(24) of the 1/f noise intensity (25). 

V. CONCLUDING REMARKS 

Analysis of the exactly solvable model of 1/f noise display main features of the noise and 
serves as a basis for revealing of the origin and parameter dependences of the flicker noise. 
This allows us to make generalizations of the model resulting in the expression for the 1/f 
noise intensity through the integral of the characteristic function of the distribution of the 
time intervals between the subsequent transition times of the elementary events, pulses or 
particles. 

It should be noticed, however, that Eq. (25) represent an idealized 1/f noise. The 
real systems have finite relaxation time and, therefore, expression of the noise intensity in 
the form (23) is valid only for / > (27rr re /)~ 1 with r re / being the relaxation time of the 
period's fluctuations. On the other hand, due to the deviations from the approximation 
tk+ q — tk = Tkq at large q, for sufficiently low frequency we can obtain the finite intensity of 
1/f 5 (5 ~ 1) noise even in the case ip (0) = but for the signals with fluctuations resulting in 
the dense concentrations of the transit times tk- Generalizations of the approach (21)-(25) 
and analysis of the deviations from the idealized 1/f noise expression (25) are subjects of 
separate investigations. 
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Captions to the figures of the paper 



FIG. 1. Power spectral density vs frequency of the current generated by Eqs. (2)-(4) with 
the Gaussian distribution of the random increments {Ek} for different parameters f , a, and 
7. The sinuous fine curves represent the averaged over five realizations results of numerical 
simulations, the heavy lines correspond to the numerical integration of Eq. (16) with g (q) 
from Eq. (13), and the thin straight lines represent the analytical spectra according to Eqs. 
(17) and (18). 

FIG. 2. Same as in Fig. 1 but for the uniform distribution of the random increments 
{sk}- Note that for the non-Gaussian distributions of the random perturbations we have no 
explicit expression analogous to Eq. (16) for the integral representation of the noise power 
spectral density. 

FIG. 3. Same as in Fig. 1 and Fig. 2 but for the asymmetric \\ distribution of the 
random increments {Ek}- 



10 




Fig. 1 




10" 5 10" 4 10" 3 10" 2 10" 1 f 10" 6 10" 5 



Fig. 2 




Fig. 3 



